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PROCEDE POUR L' ESTIMATION 
DBS CARACTERISTIQXJES D'X3NB PRECIPITATION 

La prSsente invention concerne un precede pour 
1' estimation des caracteristiques d'une precipitation, et en 
particulier du tavix precipitant pour une precipitation solide. 

On connait dans I'^tat de la technique !• utilisation des 
caractSristiques d'une image radar pour estimer les 
caracteristiques pluviometriques . En particulier, le brevet 
europeen EP1049944 d^crit une technique pour 1 ' estimation de 
pluie avec \xa radar. Selon ce brevet, on realise les Stapes 
suivantes : 

- on mesure , au moyen dudit radar bipolaire , sur un 
intervalle [ro, ri] donnS de rayon de trajet par 
rapport audit radar, la phase dif f Srentielle «&dp et 
la reflectivity apparente Z. selon au moins 1 ' une des 
polarisations H ou V ; 

-on determine une estimation de la valeur No* 
representative de la distribution dimensionnelle des 
gouttes de pluie, a partir de la difference de la 
phase differentielle entre ro et ri et partir d'une 
integrale d'une fonction de la reflectivite apparente 
Z, le long de 1' intervalle [ro, ri] ; 

- on deduit la valeur du taux de precipitation en un 
point si partir de No* et de la reflectivite apparente 
en ce point . 

On connait egalement un precede pour 1" estimation de 
pluie decrit dans le brevet PCT WO03007016 decrivant un 
precede pour 1' estimation d'un taux precipitant au moyen d'xin 
radar bipolaire, caracterise par les differentes etapes 
suivantes : 
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- on mesure au moyen dudit radar bipolaire, sur un 
intervalle [rl, rO] donn^ de rayon r de trajet par 
rapport audit radar, la phase dif f Srentielle Odp et 
la reflectivity attenuee Z selon au moins I'une des 

5 polarisations H ou V ; 

- on determine une estimation de la valeur K(rO) de 
1 ' attenuation specif ique en rO i. partir du profil de 
r^flectivite attenuee ainsi mesurg, ainsi qu*a partir 
de la difference de la phase dif f erentielle entre rO 

10 et rl; on determine une estimation K(r) de 

1 ■ attenuation specif ique en r en fonction de 
1 ' attenuation K(rO) ainsi determin^e et du profil de 
rSflectivitS attgnuSe Z(r); on determine le taux de 
precipitant R(r) connaissant K(r) . 
15 Ces differentes solutions permettent de caracteriser des 

precipitations liquides et d'estimer de fagon precise le taux 
de pluie (en mm/h) , mais pas d'estimer les caracteristiques 
des precipitations solides telles que la neige. 

L' invention concerne un nouveau precede visant ^ remedier 
20 ^ cet inconvenient. 

A cet effet, 1 ' invention concerne selon son acception la 
plus generale un procede pour 1 'estimation des 
caracteristiques d'une precipitation, et en particulier du 
taux precipitant pour une precipitation solide comportant une 
25 etape d ' acquisition d'une image radar comportant au moins un 
plan vertical d'une zone de precipitation et d'un traitement 
d'un profil vertical pour fournir des signaux numeriques 
representatif s de la reflectivite selon la direction verticale 
h caracterise en ce que I'on realise une etape d' integration 
30 desdits signaux representatif s de la reflectivite pour 
deiivrer un signal representatif du profil dans le plan 
vertical du diametre moyen des particules pondere par la masse 
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de chaque particule, et une 6tape de determination de la 
concentration des particules solides k partir des signaux 
calculus dans les Stapes prScSdentes. 

De preference, I'etape d' integration consiste k 
5 determiner la variable Z (h) 1' observable radar en mm^/m^ en 
fonction de 1' altitude h a partir de ladite image radar, et a 
determiner ledit diamdtre moyen des particules Dm(h) par 
resolution de 1' equation : 



■ Z est 1' observable radar k inverser en mm^m"^ ; 

■ D„ est en m ; 

■ a et Jb sont des coefficients specif iques des 
particules de type « agregats ». Le coefficient a est par 

15 exemple ggal a 35184 et le coefficient b est egal a 3.16. 

■ k^ff est le coefficient d'efficacite du processus 
d'agregation a ajuster, ce coefficient k^^^ etant par 
exemple Sgal k 0.3 

L' integration de (2) ngcessite une condition a la limite 
20 d' integration. Avantageusement , ladite condition k la limite 
d' integration est determinee afin que la valeur D^(h) au sommet 
du nuage corresponde a la valeur predeterminee du nombre total 
de particules au sommet du nuage. 

Selon un mode de realisation avantageux, on determine le 
25 profil du nombre total de particules n^ (h) par 1' equation 
suivante : 




(2) 



10 



nT(h)=x.Z(h) /D,(h)^ 

ou X est egal a 25,4 10"" 

Selon un autre mode de realisation, on determine le 



30 



param^tre meteorologique Ng (h) par 1' equation suivante : 
No(h)=y.Z(h)/D„(h)^ 
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Ou y est egal k 102 10"" 

Selon un troisifeme mode de mise en oeuvre, on determine 
le paramStre mStgorologique correspondant au profil du contenu 
en glace INC(h) [en g/m^] par 1' Equation suivante : 

IWC(h)=wZ (h)/D„(h)^ 

oH w est egal & 1,25 10"" 

Selon un quatri&me mode de mise en oeuvre, on determine 
le paramStre meteorologique correspondant au profil du taux 
precipitant R{h) solide {mm/h Equivalent fondu) par 1' Equation 
suivante : 

R{h)«r,Z{h) /D,(h) 

Oa r est egal a 4,698 lO"^*" 

L' invention sera mieux comprise a la lecture de la 
description qui suit, se refErant §l un exemple non limitatif 
de realisation. 

- La figure 1 represente un exemple de profil 
vertical de Z & inverser (dans cet, exemple, I'isotherme 0°C 
est au niveau du sol) . 

- La figure 2 represente la comparaison du profil 
de Dm resultant de 1' inversion de Z par le modele d'agregation 
avec I'estimateur classique. 

La figure 3 represente la comparaison des 
profils de NO et nT resultant de 1' inversion de Z par le 
modele d'agregation avec I'hypothese classique et des 
observations . 

- La figure 4 represente le profil vertical de Z a 

inverser 

- La figure 5 represente la comparaison du profil 
de Dm resultant de 1' inversion de Z par le modele d'agregation 
avec I'estimateur classique. 

- La figure 6 represente la comparaison des 
profils de NO et nT resultant de 1' inversion de Z par le 
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modele d'agrggation avec I'hypoth^se classique et des 
obseirvations . 

- La figure 7 represents la comparaison du profil 
de R resultant de 1' inversion de Z par le modele d' agr^gation 

5 avec I'estimateur classique. 

- La figure 8 reprSsente la sensibilite de la 
restitution du ta\rx: precipitant equivalent (mm/h) a la loi de 
densitS des particules p (D) D-y. 

- La figure 9 represente la sensibilite de la 
10 restitution de la vitesse Doppler verticals a la loi de 

densite des particules p(D) «= D-y* 

Le proc€d6 s' applique essentiellement aux precipitations 
stratif ormes . II considers que les noyaux glagogenes ne sont 
actives qu'a temperature fortement negative, c'est-^-dire au 
15 soramet du nuage. Les cristaux de glace formes a haute altitude 
sedimentent et grossissent au cours de leur chute soit par 
sublimation de la vapeur d'eau saturante ambiante, soit par 
collection et givrage de gouttelettes d'eau nuageuse 
surfondues, soit par agregation au hasard de leurs collisions 
20 avec d'autres cristaux de glace. Des trois processus de 
croissance, seule 1' agregation change la concentration en 
particules de glace. Le precede de caracterisat ion repose 
essentiellement sur une description simplifiee du mecanisme 
d'agr6gation. L' invention est basee sur une classe de precedes 
25 « profileurs », en ce sens qu'elle inverse le profil vertical 
de reflect ivite mesure dans la glace, pour en tirer le profil 
vertical du taux precipitant solide. 

Les Stapes de la methode d' inversion sont les suivantes : 
1- La distribution dimensionnelle des particules 
30 exprimSe en « diamfetre equivalent fondu est suppos6e 

exponentielle i.e. : 

N{D) =Noexp (-4D/Dm) (1) 
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oH N{D) est la concentration en particules par et 
par intervalle de diamStre, et No et D„ sont les deux 
paramdtres qui caractSrisent la distribution* 

2- On determine le sommet h^^ et la base h„i« de la 
couche de precipitation solide. 

-hmax est 1' altitude maximum du profil de 
reflectivity mesur6 Z{h) . 

b. est soit 1' altitude de I'isotherme 0°C si 

la temperature au sol est positive, soit le niveau du 
sol si la temperature au sol est negative. 

3- On determine alors le profil du paramStre D^{h) 
entre h^^^ et h^^^ en resolvant 1' equation 
dif f erentielle : 



od : 

Z est 1' observable radar k inverser en 
mm^m"^ ; 

D„ est en m ; 

a et i? sont des coefficients specifiques 
des particules de type « agr^gats », egaux 
respectivement a 35184 et 3.16 a partir des 
observations de Locatelli et Hobbs (1974) ; 

k^ff est le coefficient d'efficacite du 
processus d'agr^gation k ajuster (la valeur k^f^ 
= 0.3 semble correcte) . 
4- L' integration de (2) s'effectue k partir du 
haut, ou la condition k la limite s'exprime en fixant le 
nombre total de particules (ou nombre de noyaux 

glagogSnes actives au sommet de nuage) . On peut prendre 
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■Ktr (^inax)= ^"^i qu^l permet d'exprimer la condition a 

la limite D„(h,nax) cottime : 

Dm(Hma:^^25J lOr'^ (Z(hma.)/m(hma^f^ ( 3 ) 

5- Le profil D„(h) entre h„ax i^min ^^e fois 

5 determine, on calcule les profile des autres param^tres 

d'intSrSt par les expressions suivantes : 

a. Profil de : No (h) =102 . lO'^Z (h) /D„{h) 

b. Profil du nombre total de particules xi^i^h) 
[en m-^ : N^di) =102 . lO'^Z (h) /D^(h) ' 

10 c. Profil du contenu en glace JPirC(h) [en g/m^] 

: IWC(h)=l,25.10'-"Z(h) /D^(h)^ 

d. Profil du taux precipitant Rkh.) solide 
(mm/h equivalent fondu) 

En utilisant la loi de vitesse terminale de 
15 chute determinSe par Locatelli and Hobbs pour des 

agrSgats: [Vj. =107.6 (D in m) ] , s'exprime 

par : 

R (h) =4,698, 10""Z (h) /D^ (h) ^'^^ 
La description qui suit concerne une autre version du 
20 precede de traitement pour la determination du taux 
precipitant . 

Pour la pluie, 1 ' algor i thme utilise est 
I'algorithme ZPHI objet des brevets FR9800714 et FR0109206 • 

Pour les precipitations solides, le nouvel 
25 algorithme qui vient d'etre developp§ est 1' objet de la 
description qui suit. 

Get algorithme d' estimation du taux precipitant 
pour une precipitation solide fait partie de la classe des 
algorithmes « profileurs », en ce sens qu' il inverse le profil 
30 vertical de reflect ivite mesure dans la glace, pour en tirer 
le profil vertical du taux precipitant solide • 



I 
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L' algorithme s' applique essent iellement aux 
precipitations stratif ormes . II consid^re que les noyaux 
glagogenes ne sont actives qu'& temperature fortement 
negative, c'est-S-dire au sommet du nuage. Les cristaiox de 
5 glace formes & haute altitude s^dimentent et grossissent au 
cours de leur chute soit par sublimation de la vapeur d'eau 
saturante ambiante, soit par collection et givrage de 
gouttelettes d'eau nuageuse surfondues, soit par agr^gation au 
hasard de leurs collisions avec d'autres cristaux de glace. 
10 Des trois processus de croissance, seule I'agregation change 
la concentration en particules de glace. L' algorithme repose 
essentiellement sur une description simplifige du mecanisme 
d'agrSgation. Les 6tapes de la mSthode d' inversion sont les 
suivantes : 

15 La distribution dimensionnelle des particules 

exprim^e en « diamStre equivalent fondu est supposee 

exponentielle i.e. : 

NiD)==Noexp(r4D/nn) , (1) 
ou: 

• D est le diam^tre equivalent fondu de la particule 
de glace, 

• N(D) est la concentration en particules par m^ et 
par intervalle de diametre, 

• et No et sont les deux paramStres qui 
caractSrisent la distribution. 

On determine le sommet h^^,, et la base h^^^ de la 
couche de precipitation solide. 

30 

• h^ax est 1» altitude maximale du profil de 
reflectivity mesurS Z(h) . 



20 



25 



10 



25 
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• Kixn ®st soit 1' altitude de I'isotherme O^C si 
la temperature au sol est positive, soit le 
niveau du sol si la temperature au sol est 
negative • 

On determine alors le profil du paramdtre D„(h) 
entre et en resolvant 1' Equation dif f erentielle, dont 

la solution peut Stre determin§e analytiquement : 



(2) 



oil : 

J(M)= nexp(-4jc)x* r^-JC*"^ Jxp(-4jc0d5c'^ d^crit la 
Jo Jo 

frequence de collision 
15 • Z est 1 ' obsearvable radar ^ inverser en mm^m"^ ; 

• Da, est en m ; 

• a, Jb, c et d sont des coefficients dependant de la loi 
de densite F(D) des particules de glace, supposee 

20 varier en D""^ (ou y peut prendre des valeurs entre 0.25 

et 1.1 suivant le type des particules). Ces 
coefficients sont tires de la theorie de Mitchell 
(Journal of Atmospheric Sciences, 53, 12, 1996) pour 
representer les lois puissance suivantes : 



= cjD^ {oxX Vt est la vitesse terminale de 
chute de la particule de diametre Equivalent fondu D) 

ji = aD^ (ou A est la section ef f icace de la 
particule de diametre equivalent fondu D) 



30 
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Avec des coefficients dependant de la frequence et 
d'une fonction de (D) , repr^sentant par une loi puissance la 
relation entre la reflectivity (parametre mesur6) , la 
concentration des particules et leur diam^tre Equivalent 
5 moyen : 

Z = dNoDj 

kgff est le coefficient d'efficacite du processus 
d'agregation a ajuster (la valeur k^ff = 0.3 seinble correcte 
juste au-dessus de I'isotherme O^C) . 
10 L' integration de (2) s'effectue partir du haut, 

ou la condition a la limite s'exprime en fixant le nombre 
total de particules (ou nombre de noyaux gla9ogenes actives 

au sommet de nuage) . On peut prendre n^. (iimax)= 10^ m*^, ce qui 
permet d'exprimer la condition a la limite D„(h,^) comme : 



15 



Le profil D„{h) entre h^^^ et h^^^^ une fois 
20 determine, on calcule les profils des autres paramStres 
d' interet par les expressions suivantes : 

Profil de Nq : 



25 Profil du nombre total de particules n^ih) [en m'^] 

nrih) = 025Noih')D„ih) 

Profil du contenu en glace IWC{h) [en g/m^] 
IWC(h) = (1.2272)10^ Noih)D„ihf 
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Profil du taux precipitant R{h) solide (mm/h 
Equivalent fondu) : 

5 Profil vertical de la vitesse Doppler V^Ch) (en 

m/s) , donn6e par une loi puissance dependant de la densite 
_(D) : 

10 

Du fait de la dependance des coefficients vis-a- 
vis de la loi de density F(D), la restitution du taux 
precipitant depend de fagon critique du parametre y, comme en 
t6moigne la figure 7 . Cette loi de densite constitue done un 

15 paramdtre clef. 

En operant en tir vertical, on peut mesurer la 
Vitesse Doppler Vj^ili) , et la comparer au Vjyih) calcule par 
I'algorithme. Le profil de vitesse Doppler est discriminant, 
comme en tSmoigne la figure 8. Par ajustements successifs, on 

20 peut ainsi determiner le paramStre y de la loi F(D) pour 
lequel la coherence entre V^Ch) calcul^e et (h) mesuree est 
realisee. La figure 4 reprSsente un exemple de profil inverse 
par I'algorithme (dans cet exemple, I'isotherme O^C est au 
niveau du sol) . La figure 5 repr^sente les courbes resultant 

25 de 1 ' application du modele d' agr^gation et du module 
classique, permettant la comparaison du profil de resultant 
de 1' inversion de Z par le modele d'agrSgation avec 
I'estimateur classique. 

La figure 6 repr^sente la courbe des profils en 

30 fonction de 1' altitude (m) au-dessus de I'isotherme 0«*C. Elle 
conduit k la comparaison des profils de Nq et Jiy resultant de 
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1' inversion de Z par le module d'agregation avec I'hypothese 
classique et des observations. La figure 7 repr^sente la 
comparaison du profil de R resultant de 1' inversion de Z par 
le module d'agrSgation avec 1' estimateur classique. La figure 
5 8 reprSsente la sensibilite de la restitution du taux 
precipitant Equivalent (mm/h) ^ la loi de densite des 
particules p (D) oc D*"^. La figure 9 reprgsente la sensibility de 
la restitution de la vitesse Doppler verticale & la loi de 
density des particules p (D) «: b"^. 

10 
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REVENDICATIONS 



1 - Precede pour !• estimation des caract^ristiques d'line 
precipitation, et en particulier du taux precipitant pour une 
5 precipitation solide comportant une etape d' acquisition d'une 
image radar comportant au moins un plan vertical d'une zone de 
precipitation et d'un traitement d'un profil vertical pour 
fournir des signaux numSriques repr^sentatif s de la 
reflectivite selon la direction verticale z caract^rise en ce 

10 que I'on realise une etape d' integration desdits signaux 
representatif s de la reflectivite pour deiivrer un signal 
representatif du profil dans le plan vertical du diamStre 
moyen des particules pondere par la masse de chaque particule, 
et une etape de determination de la concentration des 

15 particules solides a partir des signaux calcules dans les 
etapes precedentes . 



2 - Procede selon la revendication 1, caracterise en ce 
que 1' etape d' integration consiste a determiner la variable 
20 Z{h) 1' observable radar en mm^/m^ en fonction de 1' altitude h a 
partir de ladite image radar, et a determiner ledit diamdtre 
moyen des particules Dm(h) par resolution de Inequation : 

^i^0.25^,^«D^nO-i«Z+|^X (2) 

o\X : 

25 ■ Z est 1' observable radar a inverser en 

mm^m"^ ; 

■ est en m ; 

■ a et b sont des coefficients 
specifiques des particules de type 

30 « agregats » • 
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■ kefc est le coefficient d'efficacite 
du processus d'agrSgation a ajuster. 

3 - Proc€d4 selon la revendication 2, caract6ris§ en ce 
que le coefficient keff est ggal a 0.3. 

4 - ProcedS selon la revendication 2 ou 3, caractSris^ en 
ce que le coefficient a est egal k 35184. 

5 - Proc6dS selon la revendication 2, 3 ou 4, caractSrise 
en ce que le coefficient b est Sgal a 3.16. 

6 - Precede selon I'une quelconque des revendi cat ions 2 a 
4, caracterise en ce que ladite constante d' integration est 
d^terminee afin que la valeur D^(h) au sommet du nuage 
corresponde il la valeur predeterminee du nombre total de 
particules au sommet du nuage. 

7 - Precede selon I'une quelconque des revendications 
precedentes, caracterise en ce que I'on determine le profil du 
nombre total de particules nt (h) par 1' equation suivante : 

nT(h)=x.Z(h)/D^(h)^ 

8 - Precede selon la revendication pr^cedente, 
caracterise en ce que x est ^gal IL 25,4 10"^®. 

9 - Precede selon I'une quelconque des revendications 
precedentes, caracterise en ce que I'on determine le parametre 
meteorologique No(h) par 1' equation suivante : 

No(h)=y.Z{h) /D„{h)\ 
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10 - Precede selon la revendicat ion prScSdente, 
caract#risg en ce que y est 6gal & 102 10"^°. 

11 - Precede selon I'une quelconque des revendicat ions 
5 prScSdentes, caracteris^ en ce que I'on determine le param^tre 

meteorologique correspondant au profil du contenu en glace 
INC{h) [en g/m^] par 1' Equation suivante : 
IWC(h)=wZ(h)/D„(h)\ 

10 12 - Precede selon la revendicat ion precedente, 

caractgrise en ce que w est egal & 1,25 10"". 

13 - ProcSdg selon I'une quelconque des revendicat ions 
prSc^dentes, caracterisg en ce que I'on determine le paramdtre 

15 mStSorologique correspondant au profil du taux precipitant 
R{h) solide (mm/h equivalent fondu) par 1' equation suivante : 
R(h)=r.Z(h)/D,{h)^'^^ 

14 - Precede selon la revendicat ion precedente, 
20 caracterise en ce que r est egal a 4,698 10'^°. 
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